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A STUDY OF POSITIVE ENERGY CONDITION IN BIANCHI V 
SPACETIMES VIA NOETHER SYMMETRIES 

SAJID ALI, IBRAR HUSSAIN 


Abstract. In this paper we use Noether symmetries of the geodesic Lagrangian in Bianchi 
V spacetimes to study various cosmological solutions of Einstein’s field equations. Our 
first result is the identification of the subalgebras of Noether symmetries of the equations 
of motion in such spacetimes with dimension 4, 5, 6, 7, 9 or 10 of the maximal algebra 
of Lie point symmetries of dimension 13. Secondly we give physical interpretation of new 
cosmological solutions which satisfy positive energy condition and yield critical bounds on 
the expansion coefficient a, in which the underlying non-flat spacetimes carry interesting 
physical properties. Specifically the energy density behaves in one of the following ways, (i) 
It is positive and constant for all time, (ii) It varies with time and attains a global maximum 
after some time and then asymptotically converges to zero, (iii) It increases for all time and 
attains a maximum value at the asymptotic limit t —> oo. In particular a non-flat spacetime 
is obtained that mimics the expansion in a flat FRW universe dominated by vacuum energy 
such that the expansion factor has the same form in both. However, the energy density is 
dynamical in the former. 


1. Introduction 

Finding exact solutions of the Einstein’s field equations (EFEs) is one of the important and an 
old problem in General Relativity (GR). The search for new solutions opens new avenues to 
our understanding of the universe. These equations are the result of Einstein’s revolutionary 
idea that the existence of matter induces curvature in a spacetime. This is responsible for an 
inward pull of all surrounding objects, thereby completely changing the standard perspective 
of gravitational force mu, Recent developments in astrophysics and cosmology reveal that 
our universe is mainly composed of a mysterious form of matter (dark energy and dark 
matter) which is the main cause of cosmic acceleration at a faster pace today [3j [4] . This 
astonishing observation laid down the basis of numerous phenomenological models and led 
to the departure from the standard solutions of EFEs in the thrust of new exact solutions 
possessing revived geometrical and dynamical properties. Our attempt in this direction is 
to exploit the geometrical features in homogeneous cosmological solutions of these equations 
which mainly arise from Noether symmetries. 

Suppose M is a four-dimensional smooth manifold equipped with a non-degenerate metric 

g of Lorentzian signature (+, —, —, —). The EFEs of the universe in standard gravitational 
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units c = 1 = G, are a system of partial differential equations G(k = 1), where G^ u 
is the Einstein tensor which contains all the basic geometric properties of spacetimes and T pv 
is the stress-energy tensor that describes the density and flux of energy and momentum in 
spacetimes [2]. The exact solutions of EFEs either arise from the geometrical consideration 
of spacetimes or by studying physical characteristics of matter or a combination of both. 
In cosmology, homogeneous spacetimes of dimension 1 + 3 are characterized using Bianchi 
classification of real Lie algebras of dimension three which resulted into nine classes I, II, ...,IX 
[5j. The famous cosmological solutions of EFEs known as Friedmann-Lematre-Robertson- 
Walker (FRW) metrics are both isotropic and homogeneous and are particular cases of 
Bianchi types I, V, VII and IX [T] . Therefore it is important to investigate other Bianchi 
homogeneous cosmologies with regard to their invariant geometrical properties. Indeed the 
authors [6], have done the symmetry analysis of all Bianchi spacetimes in the presence of a 
dynamical held. In [7], authors have classified Bianchi V spacetimes using symmetries of the 
curvature and Weyl tensors known as curvature collineations (CC) and Weyl collineations 
(WC), respectively. 

The concept of symmetries is closely related to the invariance of geometrical quantities 
under certain diffeomorphism on the manifold that could change the manifold however it 
keeps the underlying structure intact. This is a well known problem of equivalence of the 
two geometric objects under diffeomorphisms Hi- The invariance of geometrical quantities 
which are described by tensor fields on the manifold under certain diffeomorphisms is ensured 
if the Lie derivative vanishes along the vector held representing a how. Suppose T is a tensor 
held of any type on the manifold M. The key to check the invariance of certain geometrical 
quantities is to use the definition of Lie derivatives and require that TxT = 0, where +x is 
the Lie derivative along the vector held X [5]. In case of the metric tensor the symmetries are 
known as conformal Killing vectors (CKV) provided £xg^v = ^g^u, where T(x At ) is known 
as a conformal factor. If T pi/ ^ 0, then a CKV is said to be proper otherwise it reduces to a 
special CKV if = 0 and T M ^ 0. The homotheties (HV) and Killing vectors (KV) arise 
if T !A1 = 0 and T = 0, respectively. Furthermore the symmetries of the curvature tensor 
are obtained by replacing the tensor held with curvature tensor £x = 0. In order to 
investigate pure gravitational held of spacetimes Weyl tensor plays a significant role as it is 
conformally invariant and its symmetries are obtained by the same condition that the Lie 
derivative of the Weyl tensor C^ pa vanishes. 

Noether symmetries play an essential role in finding conservation laws of the equations of 
motion with the use of Noether theorem 0 EDI. Sometimes it is difficult to obtain exact 
solutions of the equations of motion, however the dynamics can be reduced by investigating 
the invariant properties of the system provided the problem under consideration is variational 
and there exist a Lagrangian. Noether’s theorem provides an explicit formula of a conserved 
quantity for each continuous symmetry transformation that leaves the action invariant [ 113 - 
In addition Noether symmetry analysis reduces the problem by specifying the unknown 
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functions that appear in the Lagrangian. This line of approach has been followed by several 
authors, notably in the pioneering works of S. Capozziello [fl2l IT3] . G. E. Prince mm 
and M. Tsamparlis mm- The connection between symmetries of the underlying manifold 
and that of the differential equations was discussed in [[18] • The classification of spherically 
symmetric static spacetimes via Noether symmetries is done in [19]. Besides approximate 
Lie symmetries was used to resolve the problem of energy in general relativity in [[201 [21 ] . 
In this paper our main interest is the classification of Bianchi V spacetimes using Noether 
symmetries. These symmetries provide crucial physical information about the conserved 
quantities of a physical system. 

We start our investigation by considering the geodesic motion in Bianchi type V spacetimes. 
In such spacetimes there are three arbitrary functions whose specification arise from the 
presence of Noether symmetries. We give a complete classification of Noether algebras along 
with first integrals for each case. It is found that Noether algebra of the symmetries of the 
geodesic Lagrangian in such spacetimes can have dimension, 4 — 7, 9 or 10. The algbera of 
Lie point symmetries of the equations of motion is of dimension, 5 — 8, 10, 12 or 13. The 
connection between Lie and Noether symmetries with symmetries of the spacetimes like HVs, 
projective collineations (PCs) and CCs has already been discussed in [18|[T6]. We highlight 
the important features of new cosmological solutions in the light of above results. 

It is well known that there are three cosmological models of our universe which are filled 
with vacuum energy, radiation and matter [I]. The large astrophysical data suggest a widely 
accepted view that our universe is nearly flat. It would be interesting to investigate as to 
what degree we can relax the condition of flatness of the underlying manifold such that the 
non-flat spacetime still carry all the important features of a realistic and viable cosmological 
model. It turns out that using Noether symmetries we can specify certain non-flat spacetimes 
which possess very nice physical properties. In particular we determine a non-flat spacetime 
that mimics the behavior of a flat spacetime filled with vacuum energy. Interestingly the 
expansion factor (so is the Hubble parameter) in the non-flat and flat spacetimes is same. 
The energy density is fixed in the flat model (hied with vacuum energy) and dynamical in 
the other case. Based on our finding we further investigate non-hat spacetimes and compare 
these with hat spacetimes dominated by matter (Einstein-de Sitter universe) or radiation. 

The paper is subsequently divided into three main parts. The Noether symmetry analysis 
is presented in the next Section where we give preliminaries to the symmetry approach. In 
the third Section we investigate the physical characteristics of such spacetimes and examine 
the implications of positive energy condition on the solutions. To confront our results with 
other extended theories of gravity we consider a particular Bianchi V spacetime which meets 
the positive energy condition and modify our results in /(i?)—gravity. It is found that 
f(R) oc i? 3 / 2 , for a dust cloud in an anisotropic Bianchi V spacetime. Lastly we summarize 
our results in the last Section. 
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2. Noether symmetry Analysis 


The Bianchi type V spacetimes have the form 

ds 2 = df 2 - A(t) 2 dx 2 - e ax {B(t) 2 dy 2 + C{t) 2 dz 2 ) , (1) 

where A, B, C , are three non-zero arbitrary functions of the cosmic time t. The expansion 
factor a, is a non-zero constant which has the units of inverse length. The above spacetimes 
represent a specially homogeneous and anisotropic cosmologies as the coefficients depend on 
the time variable t. Since these are non-static spacetimes and thus do not admit a time 
like Killing vector field. All Binachi V spacetimes admit a three-dimensional Lie algebra 
of spacelike KVs (corresponding to the isometry group G 3 , that acts transitively on the 
spacelike hypersurfaces) [5] 


x _ d x _ 9 
A 1 “ O - ’ A2 ~ O - ’ 

ay dz 


<9 d d 

X 3 = 2-- ay- - az- 

oy 


dx dy dz 

Subsequently we identify this algebra with /C 3 = {Xi,X 2 ,X 3 }, corresponding to the basic 
isometry algebra of the underlying Bianchi V spacetime. 


An action A = f Lds, of the geodesic motion possesses the Lagrangian L = L(s, x^, x M ), in 
the background of Bianchi V spacetimes, equivalently 

L = t 2 -A 2 x 2 -e ax (B 2 y 2 + C 2 z 2 ), (2) 


where an over dot represents the derivative with respect to the geodetic parameter s. A 
Noether symmetry is a vector fielcO 


x- .<9 o d , d 2 d 3 d 
A = + rf— + T] 1 — + rj 2 — + rj 3 4 — , 

ds dt dx dy dz 

which leaves the action invariant such that it satisfies the condition 

X^L + L{D S Z) = D S G, 


(3) 


(4) 


where £ = £(s,x M ), = p M (s, x M ), (// = 0,.., 3) and G(s,x M ) is an unknown function of the 

corresponding symmetry. The first integrals of the equations of motion are determined by 
famous Noether theorem which gives an explicit formula mi 

dJ 

I = ZL + {rf-in)—-G, (5) 

where x M = (t, x, y, z) and = (77 0 , 77 1 , rj 2 , p 3 ) denote the coefficients of the Noether symme¬ 
try ((3]). Note that D S I = 0, upon using the equations of motion, where 


D 


S 


d-d d d d 
Ts +t m +x di + y 'Fy +z Yz- 


( 6 ) 


1 Note that the action of the underlying Lie group is extended to the product manifold Mxl, which is 

(4+ 1)—dimensional so is to bring the dynamical symmetries (Lie and Noether) and geometrical symmetries 

(CKV,CC,PC etc) on equal footing. The action naturally includes the parametrization with respect to the 
geodetic parameter s, where this line of approach was followed in m- 
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The equations of motion (geodesic equations) for the geodesic Lagrangian in Bianchi V 
spacetimes comprise of a system of four second-order ordinary differential equations in the 
field variables 


i + e ax (BB' ij 2 + CC' z 2 ) + AA' x 2 = 0, 
2 A'.. ae ax (B 2 y 2 + C 2 z 2 ) 


x H-— xt — 

A 


2 A 2 


= 0 , 


2 B’ 
~B 

2 C' 
~C 


t + ax j y = 0 , 
t + ax ) z = 0 , 


( 7 ) 


where ' represents the derivative with respect to t. It is known that Noether algebra is a 
subalgebra of the Lie algebra of point symmetries of the differential equations [S], therefore 
our study yields classification of Noether algebras of the equations of motion for the geodesic 
Lagrangian in Bianchi V spacetimes. In order to find Lie point symmetries of the differential 
equations we extend the vector field (01) to the jet space of second-order where we require 
that 

X [2] E = 0 mod E = 0, ( 8 ) 


where E. is a given system of differential equations and is the second-order prolongation 
which is defined below 


X [2] 


[i] 

X + 7^ 


d 

dxv 


+ y 


[ 2 ] d 


dx 11 ’ 


( 9 ) 


[ 1 ] [ 2 ] 

where and are determined from the formulas 


[1] 

= 

[ 2 ] 

— 


dr^_ d^ 

ds ds ’ 

[i] 

Arf_ _ d^ 
ds ds 


( 10 ) 

( 11 ) 


The Lie symmetries are obtained by applying the operator ([9]) on the geodesic equations 
while replacing the second order derivative terms in the last equation from the equations of 
motion 


Since all Bianchi V spacetimes admit /C 3 therefore the minimum dimension of the Noether 
algebra of the geodesic Lagrangian can be obtained easily. Indeed it is easy to prove that 
the minimal Noether algebra is A /4 = /C 3 © {<9s}, which is of dimension four. 

Theorem 1. The minimum dimension of Lie algebra of Noether symmetries of the geodesic 
Lagrangian (H|) that leaves the action invariant in Bianchi V spacetimes is f. 

Proof. The Lie algebra of isometries, /C3, leaves the action invariant therefore all KVs are 

Noether symmetries of the Lagrangian (J2]). It is easy to see that the global action of full 
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isometry group G 3 on the Lagrangian (J 2 ]) is 

L(s, t, x, y, z ) —■* L(s, t,x + 2e 3 , e~ at3 y + e~ ae3 z + e 2 ), 

which clearly leaves the action invariant. Therefore all KVs are Noether symmetries of the 
Lagrangian ((2j). Finally Lagrangian (J2J) does not explicitly depend on s therefore invariance 
under s— translation is trivial, i.e., L(s, t, x, y, z) — > L(s + 64 , t, x, y, z). 

Theorem 2. The minimum dimension of algebra of Lie point symmetries of the geodesic 
equations & in Bianchi V spacetimes is 5. 

Proof. Since Noether algebra is a subalgebra of Lie symmetry algebra therefore A /4 is a Lie 
subalgebra. Besides the equations of motion ((7D are invariant under scaling symmetry sd s , 
therefore for all Bianchi V spacetimes the minimum dimension of Lie algebra is five which 
is identified as A /4 © {sc/,}. 

In order to identify Bianchi V spacetimes that possess minimal Noether algebra A /4 of di¬ 
mension four and other algebras we solve the determining equations for Noether symmetries 
which are obtained by splitting equation (j4j) in monomials and obtain a set of nineteen linear 
partial differential equations 

6 = 0 , 4 = 0 , 4 = 0 , 4 = 0 , 4 - 2 riot = 0 , 

Pox ~ A 2 r] lt = 0, r) 0y - e ax B 2 g 2t = 0, g 0z - C 2 e ax rj 3t = 0, 

A 2 T] ly + B 2 e ax r] 2x = 0, B 2 p 2z + C 2 r) 3y = 0, A 2 g lz + C 2 e ax rj 3x = 0 
2^o A! + A{ 2rji x — 4) = 0, 

2rj 0 B' + B(arji - 4 + 277 2y ) = 0, 

2r] 0 C’ + C(ar)i - 4 + 2 r) 3z ) = 0, 

G s = 0, Gt — 2770 s = 0, 

G x + 2 A 2 ? 7 i s = 0, G y + 2B 2 e a ‘ l rj 2s = 0, G z + 2C~e ax r] 3s = 0. 

The above system can be integrated for various forms of the arbitrary functions. We use 
CAS (Maple) to categorize all the cases using the command ‘rifsimp’. It is interesting to see 
that in most cases the function Aft) turns out to be an affine function. Therefore the above 
system of PDEs uniquely characterize the affine form of Aft) which specify the forms of other 
arbitrary functions. We denote Noether algebra with A/), where i refers to its dimensions 
and the distinction of subcases is made through another subscript. We observe that the 
following cases arise and in each case we also obtain the first integrals of the equations of 
motion (J7J) using Noether theorem (|5]l . 

Case 1. (4-dimensional algebra) 

It is identified as A/ 4 , which arise for the geodesic Lagrangian in Bianchi V spacetimes where 
evolutionary functions are arbitrary and not among those that arise in the subsequent cases. 
The four independent invariants of the equations of motion corresponding to A/ 4 , include 
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three invariants for algebra /C 3 and fourth invariant is for d s , which is the Lagrangian itself, 
or equivalently 

I 1 = B 2 e ax y, I 2 = C 2 e ax z, / 3 = 2 A 2 x ~ ae ax (B 2 yy + C 2 zz), / 4 = L. (12) 

The Lie algebra of point symmetries of the geodesic equations includes all Noether symme¬ 
tries and an additional scaling symmetry C\ = sd s , which is not a Noether symmetry. Here 
we use £, to denote a Lie symmetry which is not a Noether symmetry with respect to a 
geodesic Lagrangian. Therefore, the Lie algebra is A /4 © {£ 1 } and is of dimensions five. 

Case 2. (5-dimensional algebra) 

There are three cases of five dimensional Noether algebras which are possessed by the La¬ 
grangian in Bianchi V spacetimes with 

2a. A(t) = Cit, B(t ) = c 2 t m , C(t) = c 3 t n , m h n , 

2b. C(t) = C\B(t), B(t ) h A(t) , C(t) 7 ^ const. 

2c. A(t) = ci, B(t) = c 2 e pt , C(t ) = c 3 e 7 *, f3 h 7 , 

where q, V i are non-zero constants. In the first case both powers m and n can not be equal 
to zero or one, simultaneously. Similarly the powers (3 and 7 can not be zero, simultaneously 


in case (2c). In the first case we obtain two Noether algebras 

J^ 5 ,a = A /4 © { 2 asd s + atd t + 2(1 — m)d x + a(m — n)zd z }, (13) 

A/ya = A /4 © { 2 asd s + atd t + 2(1 — n)d x — a(m — n)yd y }, (14) 

A4,6 = A4 © {c\zdy - yd z } , (15) 

A/ 5 .C = A4 © {9f - j3ydy - ~jzd z } , (16) 

Besides invariants A — Z 4 , an additional invariant for each case is 

h,a = ~e ax (((m — n)z — sz)c\zt 2n — c\sy 2 t 2rn ) + ——-+ tt — 1+ 

s (c 2 f 2 x 2 — i 2 ) , (17) 

ha = e ax (({ m ~ n )y + sy)c\yi 2m + c\sz 2 t 2n \ -\— -- c\t 2 x + tt — 1 + 

s (c 2 f 2 i 2 + t 2 ) , (18) 

however for other two cases the invariants include 

h,b — B 2 c\e ax (yz — zy ), (19) 

/ 5 , c = t + clf3yye ax+2pt + c 2 7 ^ie Qa:+27t . (20) 


The algebra of Lie point symmetries is six in all cases where the identification is AAyffijAi = 
sc/,}, where / G {a, b, c}. 

Case 3. (6-dimensional algebra) 

We obtain two cases in which the geodesic Lagrangian admits a 6—dimensional Noether 




algebra. These are identified as 

A4,a = A4.& © {2 asd s + atd t + 2(1 — n)<9 x } (21) 

,A4,& = A4,6 © {2ci«sa s + a(cit + c 2 )<9i + 2cicfj;} , (22) 

up to a redefinition of constant in case 2, which are possessed by the Lagrangians in Bianchi 
V spacetimes 


3a. A(t) = c±t, B(t ) = c 2 t n , C(i) = c 3 t n , n ^ 1 . (23) 

3b. A(t) = cif + c 2 , -B(t) = c 3 , C(t) = c 4 , (24) 

respectively. The invariants I\ — J 5 are same as before with an extra invariant 

= 7^4- li Clt ^ + (tt- 1 + s(c\t 2 x 2 - i 2 )) + se ax t 2n (c 2 2 y 2 + c 2 z 2 ) , (25) 

a 

= ac 2 t + ci (a(s(clx 2 — t 2 ) +ti — 1 ) — 2 c 2 i) + c 2 t( 2 c 2 i + fir) (sax — 2 )+ (26) 

cias(c 3 ?/ 2 + c 2 z 2 )e ax . (27) 

We find that there are seven Lie symmetries for the geodesic equations where we get a 
decomposition of the extra vector fields obtained in each case. In particular, the vector field 
2asd s + atd t + 2(1 — n)d x , divides into two independent Lie symmetries 

Cl ' a = s d~s ’ £2 ’ a = °^ + 2(1_n) d4’ (28) 


such that there linear combination is a Noether symmetry but independently these are not 
Noether symmetries. It is easy to check that the algebra A4,& © {£i, a , £ 2 ,a}) is a seven¬ 
dimensional Lie algebra of point symmetries for equations of motion in case (3a). For the 
other case (3b) the Lie algebra of point symmetries is A4,b © {£1 £ 2 ,&}, where 

c, ' b=s iTs’ £ « =(ci ‘ +c2 4 + 4 1 ^' (29) 


Case 4. (7-dimensional algebra) 

The geodesic Lagrangian with all evolutionary functions being equal 

A(t) = B{t) = C(t) ^ cit + c 2 , 
admits A4,& algebra along with two additional symmetries 

x. = 4 - ^4 + © (4e-“ + a V - **)) A . 


/ 2 2U d a V Z d 


(30) 


dx 2 (9?/ 4a 

• Yl = V i + 4 - “' (r - ^ Ty - - 2 & ’ 

therefore A /7 = A4,6 © {A 6 , A" 7 }. The two new invariants are 

J 6 = A 2 (4 z + 4azx + a 2 e ax ((y 2 — z 2 )z — 2yzy)) , 

I 7 = A 2 (Ay + Aayx — a 2 e ax (( y 2 — z 2 )z — 2yzz)) . 

The Lie algebra of point symmetries is A/" 7 © = sd s }, which has dimension eight. 


(31) 

(32) 





Case 5. (9-dimensional algebra) 

In this case there appear three subcases in which all evolutionary functions are specified. In 
the first two subcases, A /7 is a subalgebra of A/g, however in the last subcase only A/ 5,6 is a 
subalgebra. 


5a. It is identified as A/g, a = J\Tj © {A/ 8 , Ag}, where evolutionary functions are specified by 


A(t) = ci, B{t) = c 2 , C(t) = c 3 . 


(33) 


Additional symmetries include 


\- - 9 X - A 
X s - Ft , 


where the gauge term is constant in all cases except for Xg, which is obtained with non¬ 
constant gauge G = 2 1. The corresponding conserved quantities are 


/§ — t, Ig — St — t . 


(34) 


The Lie symmetry algebra of the corresponding geodesic equations is of twelve dimension 
which is identified as A /7 © {X 8 , X 9 } © {A% }, where A% is 


£1 



Co — t 


d_ 

dt' 


C?, — t 


d_ 

ds 


(35) 


5b. The subalgebra is A/" 9] 6 = A /7 © {X 8 , X 9 }, for 

A(t) = c\t + c 2 , B(t) = c 2 A(t), C(t ) = c 3 A(t). 

Additional symmetries include 

Cit 2 + 2 c 2 t 


(36) 


Xs — S 77-h 

OS 


2 9 s(Cit + C 2 ) 9 


Cl 


dt ’ 


G = 


2ci 


x 9 = 2 e G = Cl> 

as ci at 


with invariants 


J 8 = cis 2 (cit + c 2 ) 2 (c 2 ?/ 2 + C4i 2 )e QX + c \ ts 2 x 2 {cit + 2c 2 ) + 2c 2 (st - t)+ (37) 

+ ci ((c^i 2 — t 2 )s 2 + 2stt — t 2 ) , (38) 

I g = Cis(cit + c 2 ) 2 (c 3 y 2 + c 2 z 2 ) e ax + c 2 tsx 2 (cit + 2 c 2 ) + Ci((c 2 x 2 - t 2 )s + tt ) + c 2 t . (39) 

Now in this case the algebra of Lie point symmetries is only ten dimensional unlike the 
previous case (5a), where the extra symmetry generator arise from the decomposition of Xg, 
i.e. 

The Lie algebra A /7 © {A/ 8 } © {£ 1 , £ 2 }, is closed which is of dimension ten. 
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5c. It is identified as A/g, c = A/ 5,6 0 A 4 , for 


A(t) = ci, B(t) = c 2 e p , C(t) = c 3 e p , /3 7 ^ ± 

where A% refers to four additional Noether symmetries 

d d d 
Xe= dt~ ^ y !hj ~ ^ z dz ’ 


a 

2ci 


X 7 = 


as d 


d 


2 13 dt <9x ’ 


G = x + 


at 


2 


2 ’ 


X* = 


d 


+ 4 + 4 i (( “ 2 - 4c5/?2)(c - 2 - c ^ 2)+4 *~*■>£ 

y^(a 2 — 4c?/5 2 ) 9 


2a 


<9^ ’ 


X Q = 


—I + 4 + 4<(“ 2 - 4 © 2 )(© 2 - 4 v 2 )+ 4 cy -~- 2 «)|- 

a at ax 4ac2 dz 

yz(a 2 - Ac 2 f3 2 ) d 


2a dy 

The corresponding invariants of the geodesic motion are 

Ie = i + Pe ax+m {clyy + d^i), 

I 7 = 2 /3 c 2 (sx — x) + a(st — t), 

J 8 = (4c 2 /3 2 - a 2 )(c 2 2 y 2 y + c^(2j/z - ^))e aa:+2/3t + 4c 2 (?/(ax + 2/3i) + y), 
I 9 = (4c 2 /3 2 - a 2 )(c 2 4i + c\y(2zy - yz))e ax+2pt + 4c 2 (|/(ax + 2$) + z). 


(41) 


(42) 

(43) 

(44) 

(45) 


In this case we again obtain twelve dimensional Lie algebra of the geodesic equations which 
is A/ 5,6 © {As, X 9 } © A$, where is 


C\ = s 


d 


£2 = ( x + 


at 


ds : 

£4 = S C 3 , £5 = ( X + 


d 


2 c 2 pj ds ’ 
at \ d 2 /3 


£3 = 


2cf/3 J dt a 


d 2/3 d 

dt a dx' 

at \ d 


x + 


2cf/3 J dx 


Case 6. (10-dimensional algebra) 

There are two subcases of ten dimensional algebras of Noether symmetries corresponding to 
a = 2ci/3 and a = —2ci(3, respectively. 


6a. It is identified as A/g )C © A^, with 


A(t) = ci, B(t) = c 2 e pt , C(t) = c 3 e pt , /3 = 


- 


a 

2ci 


(46) 
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where Ay is the same as in case (5c) by substituting the value a = 2ci/3. 
symmetry A^ that arises is 


d . d t — cix d d d 

X w = 2 s— + (t- cix)—- — + y— + z— , 

as at ci ox dy dz 


with corresponding invariant 


An additional 


iio = (cix + t)(ci(si; — x) + (t — st)) + c 2 (s(c 2 y 2 + c 2 i: 2 ) — c 2 2 yy — c\zz) e 2 ^( t+c ^ x )_ ( 47 ) 

In this case we obtain thirteen dimensional Lie algebra of the geodesic equations which is 
A/'s, 6 © {A 8 , X 9 } © A^ a , where A 9a is 

. d „ ( t\d d Id 

£1 = s —, £2= X+— —, £3 = TW - W-, 

as \ ci J os dt ci ox 

_ 9 l-2j3t9 _ d 1 + 2ci j3t d 

A- S £ s , — t— H 2ci/? ~dx’ Ce ~ X m~ 2cp di ' 


6b. It is identified as A/g )C © Afd with 

A{t) = ci, £(t) = c 2 e /3t , C(t) = c 3 e^, /3 = 

where Ay is the same as in case (5c) by substituting the value a = —2ci/L 
symmetry A^ that arise is 


v d . . d t — cix d 

A10 — 2 s——I- (t + Cix) — H-— 

as at ci ax 


d d 

+ t/w- + , 

ay dz 


with corresponding invariant 


(48) 

An additional 


ho = (cix — t)(ci(sx — x) — (t — st)) + c 2 (s(y 2 + i 2 ) — yy — zz) e 2 ^ 1 C1X \ (49) 


In this case we again obtain thirteen dimensional Lie algebra of the geodesic equations which 
is A£ 5,6 © {X s , X 9 } © A^ b , where A% >b is 


£1 



£4 


s hs, 


C2 = { x -^) t s ’ c 3 

r - f d 1 ~ 9 

5 dt 2ci(3 dx ’ 


d Id 
dt ci dx ’ 

d 1 — 2ci/?t <9 
6 = Z dt 


This completes the classification of Noether symmetries of the geodesic Lagrangian in Bianchi 
V spacetimes. The equations of motion (|7J) inherit a twelve dimensional Lie algebra of Lie 
point symmetries for the simplest case, in which all three scale factors are constant which 
is not the maximal algebra as one would expect for the simplest model. The reason why 
simplest Bianchi V model does not attain a maximal algebra of dimension 13, does not 
come as a surprise because the underlying Riemannian manifold is not flat but contains 
flat sections. The algebra of Lie symmetries of the equations of motion in a flat space is 
unique and corresponds to sl(n + 2 ,R) JT 6 J. In our analysis, Bianchi V spacetimes are 

not flat in general (except when a = 0) the maximum dimension of Lie algebra of Lie 
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Table 1. Dimensions of Noether and Lie algebras of point symmetries 


Cases 

Noether Algebra J\f 

d(M) 

Lie Algebra £ 

d(C) 

1 

A/4 

4 

A/4 © {£1} 

5 

2a 

Mb,a 

5 

A/5 , a © {£l} 

6 

2b 

•A/5,ft 

5 

A/5,ft © {£l} 

6 

2c 

Mb,c 

5 

A/" 5 , c © {£l} 

6 

3a 

Me ,a 

6 

A/5,ft © {£i,cm £2,a} 

7 

3b 

■A/]),ft 

6 

A/5,ft © {£l,ft, £2,6} 

7 

4 

Af 7 

7 

A/7 © {£1} 

8 

5a 

•A/9, a 

9 

A/7 © { Ag, A 9 } © A3 

12 

5b 

A/9,6 

9 

A/7 © {Ag} © {£1, £ 2 } 

10 

5c 

-A/g,c 

9 

A/ 5,6 © {Ag, Ag} © Af 

12 

6a 

•A/10, a 

10 

A/5,6 © {A 8 , A 9 } © Af, a 

13 

6b 

A/10,ft 

10 

Mb,b © {Ag, Ag} © A^ b 

13 


point symmetries is thirteen. Noether symmetry classification of the geodesic Lagrangian in 
Bianchi V spacetimes reveal that there are six Noether algebras of dimensions 4, 5, 6, 7, 9 or 
10 which are the subalgebras of the 13—dimensional Lie algebra. Thus, we have established 
the following results. 

Proposition 1. The geodesic Lagrangian (G|) of Bianchi V spacetimes can have Noether 
algebra of dimensions 4, 5, 6, 7, 9 or 10. 

Proposition 2. The algebra of Lie point symmetries of the geodesic equations of Bianchi 
V spacetimes specified by the Noether symmetries, can have dimensions 5, 6, 7, 8, 10, 12 or 

13. 

The summary of our results is given in Table 1, where d(Af) and d(£) refers to the dimensions 
of Noether and Lie algebras, respectively. 

3. Physical Interpretation of New Solutions 

The study of inhomogeneous and anisotropic cosmologies started soon after the birth of GR 
[22] and a detailed analysis of exact solutions in terms of the asymptotic of singularities 
appeared in the works [231 [21] followed by numerous attempts m Ea m Eg. For a more 
detailed survey on such cosmologies the reader is referred to [29]. The class of Bianchi V 
spacetimes contains anisotropic and homogeneous cosmologies which are crucial to investi¬ 
gate for several reasons. The universe is homogeneous at very large scales and the question of 
whether it had started with a little bit irregularity require us to analyze small perturbations 

away from the high symmetry of Friedmann models. This could help us investigate on the 
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present-day anisotropy of the microwave background radiation and irregularities in density 
and temperature at the early epochs when the radiation was emitted. Besides adiabatic 
cooling, viscous dissipation and particle creation of an anisotropic universe can be studied 
with new cosmological models. 

The information of matter content in Bianchi V spacetimes is contained in the Riemann and 
Weyl tensors which can be used to examine the effect of tidal forces due to curvature in the 
manifold along the geodesics. The former helps us to monitor the change in volume while 
later provides change in shape of the observer along geodesics. The basic requirement for a 
spacetime to be physical is that the positive energy condition is met. The positive energy 
condition requires that T 00 which corresponds to the energy density is non-negative therefore 
Too > 0. It is worth pointing out that the existence of such Noether symmetries completely 
specify the cosmological models which is in agreement with the results in [30], where point¬ 
like Noether symmetries were employed to determine feasible models in extended gravity 
quantum cosmology. Our prime interest here is to investigate the positive energy condition 
in the models specified by Noether symmetries. 

The components of Einstein tensor, Weyl and curvature tensors involving arbitrary functions 
are already given in [7]. It was also shown that the rank of a 6 x 6 curvature matrix is 3, 
4, 5 or 6 in Bianchi V spacetimes where the case of rank 3 give rise to infinite dimensional 
Lie algebra of proper CCs. On the other hand the rank of 6 x 6 Weyl matrix is 0, 4 or 6 
and such spacetimes do not admit proper WCs except for the trivial rank zero case. The Lie 
algebra of proper conformal KVs in such spacetimes is four-dimensional. We now discuss 
the physical interpretation of the cosmological solutions obtained in the last section in the 
light of positive energy condition. 

Case 1 . (4—dimensional algebra) 

The class of Bianchi V spacetimes admitting 4-dimensional algebra of Noether symmetries 
contains all three arbitrary scale factors and the corresponding invariants are given in (fT2lh 
We investigate general physical characteristics of this spacetime by considering the energy 
density 

A' B' B' C’ A' C' 3 a 2 

Ib + Ic + Ic" 4 42 ’ 

where the positive energy condition requires that p(t) >0. It is convenient to bring p(t) into 
a more useful form 

p(t) = a'b' + b'c' + a'c' — ^ —e~ 2a , (51) 

by introducing a simple change of quantities a = In A, b — In B, c = In C. The first 
three terms in the above equation does have the qualitative behavior of kinetic energies in 
Newtonian mechanics [I]. Since a', encodes the information how fast or slow the expansion 
takes place, therefore the mixed term a'b', can be regarded as the kinetic energy of the 

composite system of both a and b. Similar argument holds for the composite systems (5,c) 
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and (a,c). However it is important to note that in the above relation we do not have the 
contribution of individual kinetic energies due to a, b and c, respectively. In fact the product 
terms can be set equivalently to 

p[t) = - [a! + b' + d) — - ( a' 2 + b' 2 + c' 2 )- —e 2a , (52) 

where the first term is positive and corresponds to the kinetic energy T c , of the composite 
system a, b and c. The second term is the total kinetic energy T e , of individual systems which 
is negative. Therefore in all Bianchi V spacetimes the total kinetic energy of individual and 
composite systems is irrelevant and the quantity which is crucial is the difference T c — T e . 
For a realistic model we require that T c > T e . The energy density is negative in the case 
when T c = T e , where the spacetime is unrealistic. The last term could be regarded as the 
potential term which is an exponential function of the scale factor a. The significance of 
the first scale factor A(t), upon the others is apparent and notably Noether symmetries also 
characterized it. 


Note that the contribution of the last term is small compared to other terms as long as |H(t)| 
grows with time. However the evolution of A{t) is critical in the interval A{t) G [0,1], in 
which case a(t) is negative and the last term plays a significant contribution in decreasing the 
overall energy density of the spacetime. On the other hand the slopes of b(t) and c(t ) play a 
significant role on the energy density of spacetimes. In particular from the equation (f5Th , it 
is clear that the contribution of the product b'c ', is larger than the sum b' + c!. Therefore if 
both b and c, simultaneously accelerate or decelerate then both slopes are positive or negative 
and we expect the spacetime to be realistic. Moreover if the slopes are opposite then the 
only possibility for a physical spacetime is that when A[t) G [0,1] , in which case the last 
term contribute in a positive energy density so is to balance the effect of other terms. Now 
we consider those cases in which the evolutionary functions are specified by the existence 
of Noether symmetries and identify the critical bounds on the expansion parameters for 
realistic Bianchi V spacetimes. 


Case 2 . (5-dimensional algebra) 

In this case we obtained three subcases in which one of the subcase contains arbitrary 
functions. We consider the cases in which all functions are completely specified 

2 a. ds 2 = dt 2 — c\t 2 dx 2 — e ax {c^t 2rn dy 2 + clt 2n dz 2 ), m 7 ^ n, (53) 


which has non-zero components of the Weyl tensor therefore it is not conformally flat space- 
time. There are eleven non-zero curvature invariants where we mention only Ricci scalar 
and denote it by I 11 , given by 



4c 2 (m 2 + mn + n 2 ) 


2 c\t 2 
14 


(54) 




which is singular at t — 0. The Einstein tensor is given by 


G 




4 cf (m+mn+n)— 3a 2 
4 eft 2 

0 

0 

a(m+n— 2) 

2 1 


0 

c 2 ((l — n)(m + n) 

0 

0 


0 

m 2 ) + a 2 /4 0 

G 22 
0 


a(m+n-2) 

2£ 

0 

0 

G33 


(55) 


where 


(a 2 — Am 2 c\)c\t 2 G ^e ax (a 2 — 4n 2 c 2 )c 2 t 2 ^ m ^e c 

G 22 =-—-j G33 =- 


4c 2 ’ 4c 2 ' (56) 

For a realistic cosmological model, we impose positive energy condition on the dynamical 
energy density 

4c 2 (m + mn + n) — 3a 2 


p(t) = 


4:C\t 2 


(57) 


which is positive if 4c 2 (m + mn + n) — 3a 2 > 0, therefore we obtain a critical bound on the 
evolution factor a 


a\ < 21ci 


m + mn + n 


(58) 


Since the quantity in the square-root must be positive therefore we obtain an extra condition 
on the powers m and n, i.e., m > —n/(n + 1), n ^ —1. If n is a positive number then there 
are two possibilities that either m is positive or negative. If it is negative then the spacetime 
becomes singular at t — 0, while it is non-singular in the other case. For — 1 < n < 0, the 
evolutionary function B(t) = Cit m , has a positive power. However in the other case n < — 1, 
the evolutionary function can have a positive or negative power as before. Therefore the 
evolution of above spacetime is such that it started at an initial time t = C 4 7 ^ 0 , then the 
energy density continue to decrease and vanishes as t —> 00 . The flux across x and y surfaces 
is zero therefore the density of first two components of linear momentum is zero however the 
z— component of linear momentum density is 


, . ^ aim + n — 2) . , 

P(t) = G„3 = --- >- . (59) 

which could be positive or negative depending on the choice of a, m and n and it asymptot¬ 
ically decays as t —> 00 . 


Similarly in the other case the spacetime has the form 

2 b. ds 2 = dt 2 - c^dx 2 - e ax (<%e 2iit dy 2 + <%e 2 ~ /t dz 2 ), 0 ± 7 , 


(60) 


which also has non-zero components of the Weyl tensor therefore it is not conformally flat. 
Again there are eleven non-zero curvature invariants and the Ricci scalar is 

3a 2 - 4c 2 (/3 2 + /3 7 + 7 2 ) 


T n _ 


15 


2 c 2 


(61) 














which is non-singular and non-dynamical unlike the previous case. The matter tensor is 
given by 


GW = 


'4/37cJ—3o 2 


4 c? 


0 

0 

q(/3+7) 


0 

c ?( 7 2 + / 3 7 + / 3 2 ) 

0 

0 


0 

(a 2 —47 2 c 2 )c|e a:c " , " 2 ^ t 

U? 




0 

0 

(a 2 —4/9 2 c 2 )c|e a3: + 2 ' y * 
4c? , 


( 62 ) 


In order to obtain a bound on the coefficient ex, we consider the density 

4/^7 c\ — 3cx 2 


p{t) = 


4c? 


(63) 


which is positive if 


cx| < 2 |ci 



(64) 


Now there are two cases that either /3 and 7 are both positive or negative. Unlike the 
previous case the above spacetime has a fixed energy density for all time. On the other hand 
the momentum density is also fixed which could be positive if both (3 and 7 are positive and 
negative otherwise. 


Case 3. ( 6 -dimensional algebra) 

The spacetime is 

3 a. ds 2 = df 2 — c\t 2 dx 2 — e ax {c\t 2n dy 2 + clt 2n dz 2 ), (65) 


where the Weyl tensor vanishes therefore it is conformally flat and Petrov type O. There are 
four curvature invariants where the Ricci scalar is 


T n _ 3(cx 2 — 4n 2 cj) 


T’<■ _ 
^ 1 ,a 


The matter tensor becomes 


G^iv — 


0 


’ 4n(n+2)c^ —3a 2 
4^72 

0 ^+nc?(2-3n) 


0 

(1— n)a 
t 


0 

0 


2 c 2 t 2 


0 

(o 2 — 4n 2 c?)c|i 2 ( n_1 U° 

u? 

0 


( 1 — n)a 
t 

0 

0 


( 66 ) 


(ck 2 — 4n, 2 c?)c 2 t 2 (" 1 ^e ax 


4c? 


. (67) 


therefore the energy density is positive G 00 = pit ) > 0 , for the following critical bound on cx 


( 68 ) 

where the above spacetime is realistic. The quantity n{n + 2) must be non-negative which 
is true if n > 0 or n < —2. The energy density vanishes as t —> oo, so is the momentum 
density for the above spacetime. 
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In the other subcase the spacetime is 

3b. ds 2 = d t 2 — ( C]t + c 2 ) 2 do; 2 — e ax (c\&y 2 + c 2 d^ 2 ), 


(69) 


where the Weyl tensor vanishes therefore it is conformally flat and Petrov type O. The 
Einstein tensor becomes 




-3 a z 

4(cit+c 2 ) 2 

0 

0 

aci 


Clt+C2 


0 

o?_ 

4 

0 

0 


0 

0 

c 2 a 2 e ax 

4(cit+C2) 2 

0 


ac\ 

C\t+C2 

0 

0 


4(cit+C2) 2 i 


(70) 


therefore the energy density is negative Goo — p < 0, for all time and the above spacetime 
is unrealistic. 


Case 4 . (7-dimensional algebra) 

In this case the spacetime involve one arbitrary function 

ds 2 = dt 2 - A(t ) 2 (da; 2 - e ax (c\dy 2 + cjdz 2 )) , (71) 

which turns out to be the case where Weyl tensor is zero thus the above spacetime is con¬ 
formally flat and is of the Petrov type O. This is an interesting case as it can be regarded 
as an inhomeogeneuos extension of the FRW spactime. There are four curvature invariants 
including the basic Ricci scalar 

3(4 AA!' - 4 A ' 2 + a 2 ) 2 


Ii = 


I n - 

J 3 — 


3 (a 2 - AAA" - A A 12 ) 
2I2 ’ 

3(47L4" - AA ' 2 + a 2 ) 3 


5127L 6 

The Einstein tensor assumes the form 


tK _ 
1 2 ~ 


T n - 

J 4 — 


64A 4 

21{AAA" - AA ' 2 + a 2 ) 4 
163844® 


(72) 


G, w = 


■3(4A ,2 -a 2 ) 

4A 2 

0 

0 

0 


0 

(a 2 -SAA"-AA' 2 ) 
4 

0 

0 


0 

0 

c\{ot 2 -SAA"-AA' 2 )e° 

i 

0 


0 

0 

0 

c%(a 2 —8AA"—4A' 2 )e ax 

i 


(73) 


where there is no non-zero component of the linear momentum density. The energy density 

3(4W 2 - ck 2 ) 


p(t) = 


AA 2 


(74) 


which could be positive or negative, in general. We now consider two dynamical evolutions 
of this non-flat spacetime in terms of a power-law and an exponential law. I 11 case of a 
power-law A{t) = cR” 1 , the density p{t) —> 00 for m < 0 and p(t ) —> 0 for m > 0, as 
t —> 00 , where the former is not an interesting case while the later contains some physical 
information. For example the density takes the form 

p(t) = 3 m 2 t - 2 - ^rt~ 2m , (75) 
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(76) 


and its behavior could be examined by choosing an initial point 

1/(1 ~ m ) 


tr, — 


{ 2mc\ \ 

V « / 


m/1 


where it is zero and it was negative before. Therefore we obtain three different cases (i) 
0 < m < 1 (ii) m = 1 and (iii) m > 1. 


In the first case we choose m = 1/2, without loss of generality and obtain 


P(t) 





tr 2 


(77) 


which is zero at t = c\/a 2 > 0. The critical point of pit) is t c = 2c\/a 2 ) where the density 
attains a global minimum because p"it c ) > 0. Since p(t c ) < 0, is negative therefore it is 
not a realistic spacetime. For the second case m — 1, the behavior of energy density is such 
that it starts from an initial value and decreases continuously till it vanishes and the global 
maximum is the initial point where the energy density has started at some non-zero time. 
The graph of energy density is depicted in Figure 1. In the last case m > 1, it attains a 
global maximum value at 


/max 


~o^~ ) 


m/1 


(78) 


after which it continue to decrease till it vanishes as t —> oo, as is shown in Figure 2. By 
introducing a Hubble type parameter responsible for the possible expansion or contraction 
of this spacetime we find that 


m = 


A' 

~A 


m 

T’ 


(79) 


therefore H(t) oc £ _1 , and this model mimics the behavior of a flat matter dominated universe 
(Einstein-de Sitter universe). Note that in which the expansion function carry the same 
form for m — 2/3. in the above case. It is interesting that our model is not flat and carries 
an energy density which is different from the energy density in a flat Einstein-de Sitter 
universe which varies as p m (t) oc f~ 2 . In our non-flat spacetime the energy density varies 
p(t) oc (4/3 — 3a 2 t 2 / 3 /4c 2 )f~ 2 in the case m = 2/3, where the convergence of p(t) is faster in 
comparison to p m . In the case of radiation dominated universe the energy density also varies 
in the same way. However it is constant for a flat spacetime dominated only by vacuum 
energy. Note that in the case of radiation dominated flat universe the expansion factor is 
proportional to t 1 / 2 , that lies in our first case. The energy density varies in the same way 
p 7 oc f -2 , which in the above non-flat spacetime is p[t) oc (1 — a 2 t/c\)t~ 2 . 


In the case of an exponential form A(t) = cie**, the energy density becomes 


P(t) 


3/3 2 a 2 2fit 

4 4c 2 
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Figure 1. The graph of Figure 2. The graph of 

energy density for a = energy density for a = 

1, m — 1, Ci = 1 1, m — 2, ci = 1 

which has the positive initial value p(0) > 0, if a < y/3ci/3. The case in which (3 < 0, 
represents a contracting non-flat spacetime that carries a positive energy for a short time 
and continue to decrease such that p —>• — oo as t —> oo, because the contribution of the 
dynamical term will dominate the first term at some point. On the other hand an expanding 
non-flat spacetime requires that (3 > 0, in which case we observe that the energy density 
increases rapidly for a short time where after it remains constant such that at an asymptotic 
limit t —» oo, it attains a finite positive value 3/3 2 /4. The graph of energy density for this case 
is given in Figure 3. It is unexpected because in an expanding flat spacetime we expect that 
the energy density decreases indefinitely, however the curvature in this spacetime confine 
it to attain a positive definite value at an asymptotic limit. The Hubble parameter for 
this spacetime is non-dynamical H(t) = (3 > 0, therefore this non-flat spacetime mimics a 
flat universe dominated by vacuum energy which has fixed density. However this non-flat 
spacetime has a varying energy density such that it attains a fixed positive value after a 
short interval of time. 

Case 5 . (9-dimensional algebra) 

Following subcases arise in which we obtain nine-dimensional Noether algebra. However the 
form of Einstein tensor can easily be obtained from the last case by substituting the value 
of arbitrary function A(t), in (1731) . 

5 (a) ds 2 = df 2 — c 2 dx 2 — e ax (c\dy 2 + c^dz 2 ), 

5 (b) ds 2 = dt 2 — (cpt + c 2 ) 2 (da; 2 + e ax (dy 2 + d^; 2 )), 

5 (c) ds 2 = dt 2 - c 2 dx 2 - e ax+pt {c 2 2 dy 2 + c 2 d^ 2 ). 
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Figure 3. The graph of 
energy density for a = 
l,(3 = l, Cl = l 


The curvature rank for the spacetimes in case (5b) is 3 and the Lie algebra of both curvature 
and Weyl collineations is infinite dimensional [7], In the first case we obtain following non¬ 
zero components of the Einstein tensor 


G 


liv 



0 0 
0 0 


0 

0 

a 2 c|e' 
4 c'f 

0 


0 

0 

0 


(81) 


Since the energy density is negative therefore the above spacetime is unrealistic although 
it is conformally flat. The components of the matter tensor in case (5b) can be obtained 
directly by placing the value of A{t) = c\t + C 2 , in case 4, equivalently 


G 


/ns 


~ 3(4 c\-a 2 ) 
2(cit+C2) 2 

0 

0 

0 


0 

a 2 —4c 2 
4 

0 

0 


0 

0 

(a 2 —4cf)e ax 
4 

0 


0 

0 

0 


(a 2 —4cf)e ax 

4 


(82) 


The positive energy condition requires that |ct| < 2|ci|, which indicates that the expansion 
coefficient a must be smaller than the twice the slope of graph of A(t) = cqf+ C 2 . The energy 
density vanishes as t —» oo. 


Lastly the case (5c) has a distinct difference from the previous two cases that there arise 
non-zero off-diagonal terms in the matter tensor. Following are the non-zero components of 
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the Einstein tensor 


G 


flU 


(3a 2 -4/3 2 cR 

4 

0 

0 

—a/3 


0 

a 2 —12/3 2 c 2 
4 

0 

0 


0 

0 


(a 2 —4/3 2 cf)c2 
4c'f 


gCx.x-\-2/3t 


0 


—a/3 


0 

0 

(a 2 —4/3 2 c 2 )c 2 e ax+2pt 
4c^ 


(83) 


where we must have |a| < 2/\/3 |ci/3|, for a realistic model. Note that the energy density 
in this Bianchi V spacetime is fixed and it attains non-zero components of the momentum 
density unlike the previous cases where these were zero. 


Case 6. (10-dimensional algebra) 

There are two subcases of 10-dimensional Noether algebra which are 


6(a) ds 2 = d t 2 - cldx 2 - e 2 ^ ClX+t \c 2 2 dy 2 + c 2 3 dz 2 ), 
6(b) ds 2 = dt 2 - c 2 dx 2 - e 2 ^~ cix+t \c 2 2 dy 2 + c 2 dz 2 ). 


Although the above models are conformally flat but these are unrealistic spacetimes as the 
Einstein tensor is 


' —2/3 2 0 0 ^2ci/3 2 " 

0 -2 cj/3 2 0 0 

0 0 0 0 

_=F2ci/3 2 0 0 0 


(84) 


and the energy density in both sub cases is negative, G 0 o = —2 f3 2 < 0. 


4. Some Considerations into /(i?)—G ravity Theory 


We now develop a few considerations for the anisotropic cosmological solutions that are ob¬ 
tained as a consequence of the existence of Noether symmetries in the context of /(i?)—gravity. 
The idea is to probe into earlier epochs in the evolution of the Universe using the results 
obtained in the previous section. As there arises several cases therefore we confine ourself to 
a few most relevant cases for our purpose. For example let us consider the case (5b) which 
meets the positive energy condition for (|a| < 2|ci|) given by 

ds 2 = dt 2 - ( Cl t + c 2 ) 2 {dx 2 + e ax (dy 2 + dz 2 )), 

in which case the Ricci scalar is 

3(« 2 - 4c 2 ) 

1 2 (at + c 2 ) 2 ’ 

which remains negative for all time due to positive energy condition therefore the underlying 

manifold is hyperbolic. We now obtain the field equations using the action of f(R)— gravity 
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involving a matter term in standard gravitational units [3TJ [32] 

A = / d 4 x^/=g f(R) + S M (g^,fp), (87) 

where k — 1, g denotes the background metric and f(R) is a general function of first 
curvature invariant which in our case is given by ()86lh Here ijj denotes all matter fields 
coupled to gravity. The held equations now assume the form 

f'R, v - l -f{R)g,v - (V,V, - 9ta , □) f = T, v , (88) 

where f = df(R)/dR , is the covariant derivative with respect to the Levi-Civita connec¬ 
tion of the metric and □ = V^V At . An immediate consequence of the above held equations 
is that their trace yield an important equation 

Rf ~~ 2/ (R) + 3 □(/ , ) = T, (89) 

which indicates how matter part T, in /(i?)-gravity is differentially connected to the curva¬ 
ture in the given spacetime unlike in standard GR, where we have an algebraic relationship 
R = — kT. We now employ Bianchi V spacetime d85ll and solve the helds equations in a 
relatively simple case of a dust cloud. As in our case the Ricci scalar R is a function of t, 
therefore the held equations yield a system of ordinary differential equations f'(R(t )) = 
where we subsequently represent the time derivative with an overdot (not to be confused 
with an overdot in Section 1). In this case the non-zero components of Ricci tensor are 

4c? — o? 

Rn =-—-, R 22 = R 33 = e° x R.\i , (90) 

which upon using the held equations (j87l) result into 

- \f(K) = Ho. (91) 

(4c? - a 2 )T - tr,t + c 2 ) 2 ( 2T - f(R)) = T„. (92) 

We now assume a perfect fluid whose energy-momentum tensor is given by 


T fl „ = (p + - pg^, (93) 

where p{t) and p(t) are energy and pressure densities of the fluid which satisfy the equation 
of state 

p = up, 0 < u < 1, (94) 

while u n = yd 7 oo(l, 0, 0, 0), is the four velocity in comoving coordinates. The conservation 
equation yields 


, . . A B C . 

P+(P + P)|-7 + -H + 75 = 0. 


(95) 


which for the underlying spacetime 
the form 


P 


ABC 

and in case of a pressure-less fluid u = 0, assumes 

3cip 


C\t + C2 

22 


= 0, 


(96) 




whose solution involving a constant of integration C\ is 

p(t) =' (97) 

By substituting the above value in (19TT) and (192]) . we obtain an ordinary differential equation 
Tit), given by 

C i 


14c? - a 2 ) 

T ~ 2 (dt + ca ) 2 * 77 + (cR + c 2 ) 3 


= 0, 


which has a solution 


J~(t) — C 2 (cR + c 2 )^ + + C^cR + c 2 )^ — 


2C'i 


(98) 


(99) 


cR(cR + c 2 ) 

where 6*2 and C 3 are constants of integration. The exponents A + and A_, are given by 

^ ci + \/9cf — 2a 2 ci - ^/9cf-2a 2 

A+ ~ 2ci ’ A -“ 2ci ’ (1 ° 0) 

which are both real for positive energy condition. In order to find the explicit form of f(R), 
we employ the equation ([89]) to obtain 

2Ci 


^ (ci t + c 2 ) 3 ’ 

which in terms of R , using ( 186 ]) is given by 

W) « fl 3/2 . 


( 101 ) 


( 102 ) 


This is quite unexpected because the above relationship is found for /(i?)—gravity in the 
background of FRW universe and it is observed that all other models other than this remain 
obscure [3T]. Therefore, an analysis of a dust cloud in the background of an anisotropic 
and homogeneous Bianchi V spacetime ([85]) which does not correspond to an FRW universe, 
reveal the same form of f(R). Here we have demonstrated a successful implementation of 
our results obtained in the previous section in /(i?)—gravity. It would be interesting to 
investigate more physical models with different equations of state in the background of other 
Bianchi V spacetimes following the same lines. A detailed analysis of rest of the cases of such 
spacetimes and a brief comparison with other extended theories of gravity will be discussed 
elsewhere. 


5. Summary 

We used Noether symmetries to study Bianchi V cosmologies. It is identified that the algebra 
of Lie point symmetries of the geodesic equations in such spacetimes can have dimensions 
5 — 8 , 10,12 or 13. On the other hand the dimensions of Noether algebras is 4 — 7, 9 or 10. 
The presence of Noether symmetries helps us to obtain first integrals or constant of motion 
using Noether theorem corresponding to the geodesic Lagrangian. A key advantage of this 

approach is the determination of unknown functions in the Lagrangian which are specified 
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while solving the determining equations for Noether symmetries that reduces the dynamics 
significantly. In all cases 2 — 6, we explicitly specified the unknown evolutionary functions in 
these spacetimes. Besides one can use first integrals to obtain closed form of exact solutions 
of the equations of motion, i.e., geodesics in non-flat spacetimes. 

In order to interpret our results physically we used positive energy condition which holds 
for all realistic cosmological models. Interestingly it imposes constraints on the solutions 
and provide us critical bounds on the constants for which a Bianchi V spacetime is physical. 
The expansion coefficient a, in all Bianchi V spacetimes is crucial. By imposing the positive 
energy condition on the obtained solutions it turns out that we can investigate different 
possible critical bounds on a, in which the non-flat spacetimes carry interesting physical 
characteristics. The positive energy condition further imposes constraints on the spacetimes 
and we obtain cases in which the energy density behaves in one of the following ways. It is 
positive and constant for all time. It varies with time and attains a global maximum after 
sometime where after it asymptotically converges to a relatively smaller but positive value. 
It increases for all time and attains a maximum value at the asymptotic limit t —> oo. A 
brief comparison of flat models (vacuum, radiation, matter) with the non-flat spacetimes is 
given in detail. 

Another important consequence of this study is the identification of f(R) oc i? 3//2 , for a dust 
cloud in an anisotropic and homogeneous spacetime (185|) which is proven to exist in FRW 
universe using standard /(A)—gravity approach. It indicates us to briefly probe into earlier 
epochs of our Universe when it was least isotropic. Therefore those Bianchi V spacetimes 
where the positive energy condition successfully holds are the best candidate for further 
investigation with more physical models with different equations of state. 
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